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ABSTRACT 


Undersea  sensors  often  inelude  an  aeoustie  modem  buoyed  above  a  seabed 
mooring  or  suspended  beneath  a  surfaee  buoy.  In  both  oases,  a  vertioal  oable  is  subjeoted 
to  horizontal  water  ourrents.  This  thesis  examines  the  two  oases,  the  first  oharaoterized  by 
a  oable  moored  to  the  bottom  of  the  ooean  with  a  buoyant  end,  and  the  seoond  being  a 
oable  suspended  or  towed  from  a  surfaee  buoy  or  Unmanned  Surfaee  Vehiole  (USV)  with 
a  weighted  end.  The  equations  of  motion  are  similar,  as  both  oases  have  an  objeot  affixed 
to  the  free  end  of  the  oable  and  the  other  end  fixed  to  a  stable  point.  A  physios-based 
algorithm  in  MATLAB  models  the  effeots  of  drag  and  buoyanoy  on  the  oable  and 
prediots  the  steady-state  shape  of  the  oable. 
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I.  INTRODUCTION 


Autonomous  ocean  sensors  serve  a  variety  of  applieations,  ranging  from 
environmental  monitoring  to  the  deteetion  of  submarines.  Depending  on  the  purpose  of 
the  sensor,  it  may  require  an  aeoustie  modem  for  sending  data  and  reeeiving  eommands. 
A  typieal  rigging  elevates  the  modem  above  a  seabed  mooring  or  suspends  it  beneath  a 
surfaee  platform,  as  illustrated  in  Figures  1  and  2.  In  both  oases,  a  oabled  system  is 
subjected  to  water  ourrents,  whioh  can  change  the  position  of  the  oable  and  the  attached 
object.  The  drag  foroe  may  shift  the  modem  out  of  a  desired  depth,  or  may  even  be  strong 
enough  to  break  the  cable.  For  a  oabled  system,  the  effects  of  flow  forces  must  be 
understood  to  prediot  cable  behavior. 

The  purpose  of  this  study  is  to  estimate  the  steady-state  shape  of  oabled  modem 
systems  in  a  oonstant  and  uniform  horizontal  flow  field.  A  finite-element  algorithm 
models  the  effeots  of  foroes  aoting  on  the  oable,  and  prediots  the  expeoted  position  of  the 
tethered  objeot.  Given  a  description  of  the  meohanioal  system  and  flow  velooity,  a 
MATLAB  program  produoes  a  two-dimensional  plot  of  the  expeoted  oable  position 
resulting  from  the  effeots  of  cable  lift  and  drag. 

Two  ease  studies  are  examined.  The  first  is  a  oable  moored  to  the  seafloor  with  a 
buoyant  objeot  at  the  free  end  (Figure  1).  The  seoond  is  a  oable  suspended  from  a 
stationary  sea-surfaoe  buoy  or  Unmanned  Surfaee  Vehiole  (USV)  weighted  by  a  dense 
objeot  at  the  free  end  (Figure  2).  The  equations  of  statio  equilibrium  are  similar  for  both 
oases.  Each  involves  a  cable  of  some  length,  or  soope,  and  an  object  attached  to  one  end 
with  the  other  end  fixed  to  a  stable  point. 
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Figure  1 .  Simple  moored  eable  subjeeted  to  a  uniform  eurrent  profile 


Figure  2.  Simple  suspended  cable  subjected  to  a  uniform  current  profde 
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II.  BACKGROUND 


The  ocean  presents  a  challenging  environment  in  which  deployed  system 
components  are  “acted  upon  by  different  forces,  including  gravity,  tensions  from  the 
mooring  lines,  hydrostatic  forces,  and  the  steady  and  unsteady  hydrodynamic  forces”  [1]. 
Examples  of  steady  forces  include  the  water  current,  or  in  the  case  of  a  towed  array,  the 
ship’s  speed.  Unsteady  forces  may  include  cyclical  waves  where  the  force  is  constantly 
changing  direction.  The  steady  and  unsteady  forces  acting  on  the  cable  determine  tension, 
bend  angle,  and  cable  shape.  This  thesis  develops  an  algorithm  to  estimate  the  cable 
shape  resulting  from  steady  forces  imparted  by  a  constant  current. 

The  effect  of  flow  forces  on  a  cable  has  long  been  studied  by  engineers.  Previous 
efforts  have  relied  on  “tables  and  figures”  to  present  “hydrodynamic  characteristics  of 
moored  array  components  such  as  buoys,  cables,  fairings,  and  sensors”  [1].  A  classic 
example  is  Pode’s  Tables  for  Computing  the  Equilibrium  Configuration  of  a  Flexible 
Cable  in  a  Uniform  Stream  [2].  Improved  estimates  are  derived  from  computer  programs 
that  predict  the  shape  of  a  cable  under  various  conditions.  For  example,  a  FORTRAN 
program  created  by  Carlson  and  Smith  has  been  used  to  determine  the  shape  and  tensions 
of  towed  arrays  [3]. 

With  their  theory  and  methods  as  a  basis,  this  thesis  develops  a  MATFAB 
program  for  predicting  the  steady-state  shape  of  a  moored  or  suspended  cable  subjected 
to  constant  and  uniform  horizontal  flow. 
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III.  CABLE  POSITION  ESTIMATION  THEORY 


This  thesis  examines  the  two  cases  of  a  moored  cable  system  and  a  suspended 
cable  system.  Displacement  of  the  cable  from  vertical  is  a  result  of  drag  forces  imposed 
by  the  horizontal  flow.  The  equations  to  calculate  the  displacement  are  the  same  for  both 
configurations.  Both  cases  are  characterized  by  a  free  and  fixed  boundary  condition.  The 
fixed  boundary  condition  at  the  end  of  the  cable  is  the  connection  to  the  mooring  or 
platform.  The  initial  forces  on  the  free  end  are  caused  by  an  attached  ballast  or  float.  The 
moored  cable  has  a  free-moving  sub-surface  buoy  positioned  at  a  water  depth  less  than 
the  fixed  end.  The  suspended  cable  has  the  end  of  the  array  free  at  a  water  depth  greater 
than  the  fixed  end. 

A,  FINITE-ELEMENT  COMPUTATION  METHOD 

The  finite-element  method  is  used  to  estimate  the  displacement  of  a  cable  subject 
to  constant  flow  forces.  Beginning  at  the  free  end  of  the  cable,  the  method  takes  the 
forces  on  the  first  element  of  the  cable  and  accumulates  the  effective  forces  along  the 
length  of  the  cable  to  get  the  total  displacement.  With  a  constant  diameter,  the  uniform 
cable  can  be  broken  into  discrete  cylindrical  elements  of  equal  length.  The  acoustic 
modem  is  treated  like  an  additional  cylindrical  element  of  the  cable,  with  differing 
dimensions  (diameter,  length,  density). 
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Figure  3.  Drawing  of  towed  eable  system  (After  [3]) 

The  formulation  by  Carlson  and  Smith,  to  determine  the  shape  of  a  eable  towed 
by  a  ship,  is  applieable  to  suspended  and  moored  eables  [3].  Figure  3  depicts  a  towed 
cable  in  a  constant  flow,  with  ship’s  velocity  being  equal  and  opposite  to  flow  speed 
across  the  cable.  In  the  case  of  the  suspended  cable,  the  ship  (buoy)  is  stationary  and  the 
“towing  effect”  is  caused  as  current  flows  from  right  to  left.  For  the  moored  case,  the 
mooring  is  stationary  at  the  bottom  of  the  ocean  and  the  “towing  effect”  is  likewise 
caused  by  the  flow  of  current. 

At  the  free  end  of  the  cable,  a  ballast  or  float  maintains  a  rigid  cylindrical  modem 
at  a  prescribed  depth.  The  initial  vertical  force  vector  on  the  free  end  of  the  cable  is  the 
gravitational  or  buoyant  force  caused  by  the  displacement  of  the  volume  of  the  ballast  or 
float  in  seawater.  The  initial  horizontal  force  vector  on  the  free  end  of  the  cable  is  the 
drag  force  imparted  by  the  current  on  the  ballast  shape.  The  initial  tension  on  the  cable  is 
the  sum  of  the  horizontal  and  vertical  force  vectors.  Incidentally,  the  initial  displacement 
angle  is  that  of  the  tension  vector. 
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The  method  now  computes  the  lift,  drag,  and  buoyancy  forces  for  each  cylindrical 
element  of  the  cable,  including  the  element  associated  with  the  modem  object.  For  this 
study,  the  properties  of  the  object  are  assumed  to  be  those  of  a  typical  acoustic  modem  as 
given  in  Figure  4. 


P 


Length  27.9  cm  (11  in) 


Diameter  10.2  cm  (4  in) 


Mass  3  kg 
Voiume  2300  cm* 

Average  Density  1.32  g/  cm* 


Figure  4.  Cylindrical  modem  dimensions 

For  the  cable  and  modem  cylinder,  the  equations  to  find  the  tension,  angle, 
vertical  and  horizontal  displacement  per  unit  length  are  the  same.  Equations  1-6  from 
Carlson  and  Smith  [3]  apply  to  a  cylindrical  shape  (the  shape  of  cable  segments  and 
modem).  Equation  1  gives  the  tension  as  the  sum  of  the  vertical  force  vector  (buoyancy 
and  lift),  the  horizontal  force  (drag),  and  the  cumulative  vector.  Equation  2  gives  the 
change  in  angle,  relative  to  the  horizontal.  The  change  in  angle  is  then  added  to  the 
previous  angle  (used  in  Equation  1)  to  find  the  tilt  angle  on  the  subsequent  segment. 
Equations  3  and  4  give  the  horizontal  and  vertical  components  of  tension  associated  with 
the  tilt  angle  in  the  subsequent  segment.  A  summary  of  forces  acting  on  a  cable  section  is 
shown  in  Eigure  5.  The  effective  drag  and  lift  coefficients  in  Equations  5  and  6  are  found 
by  utilizing  the  cross-flow  principle  outlined  in  Hoerner  [4],  with  modifications  made  by 
Carlson  and  Smith  [3].  Equations  5  and  6  take  into  account  the  drag  coefficient  Co  and 
lift  coefficient  Cl  and  the  tilt  angle  relative  to  the  horizontal  axis.  Carlson  and  Smith 
prescribe  the  values  of  Co  =  1.1  and  Cl  =  0.012  from  “Measurements  with  towed  arrays 
and  cables”  [3].  These  values  are  assumed  for  the  cable  and  modem  cylinders  unless 
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otherwise  noted.  The  drag  and  lift  forees  are  then  aeeumulated  over  the  length  of  the 
uniform  array  to  the  fixed  boundary  eondition.  Table  1  lists  the  symbol  definitions  used 
in  the  equations. 


dr 

ds 

ds 


=  gA{p-  p,^)  sin  0  +  O.5pfyV^D{C^cos0-CfSind) 

=  ~  0.5pf^V^D{Cj^smd  +  CfCosd)] 


T 


C,(^)  = 

CoiO) 


dx 

-  =  COSO' 

ds 


dy  ■  n 
—  =  smo' 
ds 


sin^  ^  cos  ^  sign  (^) 


=  sin^  6'sign(6')  +  Q 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 
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Table  1.  Cable  equations  symbol  definitions 


Symbol 

Equation  Definitions  and  Units 

g 

Gravity  aeeeleration  constant  (9.8  m/s  ) 

A 

Cross-sectional  area  of  the  cable  (m  ) 

P 

Cable  density  (kg/m  ) 

Pw 

Seawater  density  (kg/m  ) 

e 

Angle  between  the  horizontal  and  the  local  tangent  to  the  cable  (deg) 

(Positive  X  axis  is  0  deg) 

V 

Flow  velocity  (kt) 

D 

Cable  diameter  (m) 

T 

Cable  tension  (N) 

S 

Distance  along  the  cable  (measured  from  the  free  end)  (m) 

ds 

Cable  step  size  (m) 

Co 

Drag  coefficient  of  a  circular  cylinder  normal  to  flow 

Cd 

Drag  coefficient  of  an  inclined  circular  cylinder 

Fd 

Drag  force  (N) 

Cl 

Lift  coefficient  of  a  circular  cylinder  normal  to  flow 

Ci 

Lift  coefficient  of  an  inclined  circular  cylinder 

X 

Horizontal  displacement  (m) 

y 

Vertical  displacement  (m) 

signf(9; 

Signum  Function  (if  6  positive,  sign(0)  =  1 .  If  0  negative,  sign(0)  =  -  1 . 

If0  =  O,  sign(0)  =  O.) 
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Figure  5.  Forces  acting  on  a  cable  segment  (After  [3]) 

At  the  free  end  of  the  cable,  the  float  or  ballast  is  composed  of  different  materials 
depending  on  the  application.  If  the  free  end  is  designed  to  suspend  the  modem  at  a  depth 
below  the  surface  platform,  the  ballast  material  is  assumed  to  be  corrosion-resistant 
Nitronic  SOB  stainless  steel  with  material  properties  from  Bauccio  [5].  The  density  of 

-5 

Nitronic  SOB  is  calculated  to  be  7807  kg/m  by  adding  up  the  individual  densities  from 
Shackelford  of  the  component  materials  of  Nitronic  SOB  (with  the  exception  of  nitogen) 
[6].  If  the  free  end  is  designed  to  elevate  the  modem  at  a  depth  above  the  fixed  mooring, 
a  float  made  of  BZ-26  syntactic  foam  is  assumed.  BZ-26  foam  is  made  by  Engineered 
Syntactic  Systems  and  has  properties  given  in  Appendix  B  [7].  Syntactic  foam  retains  its 
buoyancy  properties  by  having  minimal  volume  change  as  pressure  (depth)  increases.  At 
a  service  depth  of  2,200  m,  BZ-26  has  a  maximum  density  gain  of  3%  over  a  24-hour 
period  [7]. 

The  algorithm  developed  in  this  thesis  allows  for  a  spherical  body,  streamline 
body,  or  no  body  at  the  free  end  of  the  cable.  The  dimensions  and  material  properties  of 
the  body  are  inputted  to  obtain  the  vertical  and  horizontal  forces  expected  on  the  shape  at 
a  particular  flow  speed.  The  methods  for  finding  the  forces  caused  by  each  shape  are 
similar,  with  differences  arising  from  the  drag  coefficient  Cd  and  cross-sectional  area  A. 
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The  streamline  body  displays  much  less  drag  than  a  sphere  for  a  given  flow,  fluid  and 
volume.  The  sphere  and  streamline  body  shapes  were  chosen  as  these  are  typical  shapes 
used  as  ballasts  and  floats.  If  no  body  is  selected,  the  initial  horizontal  and  vertical  forces 
at  the  end  of  the  cable  can  be  manually  inputted  into  the  program.  At  the  conclusion  of 
the  program,  a  two-dimensional  (2-D)  plot  in  the  x-y  plane  renders  the  shape  of  the  array. 
The  tension  and  tilt  angle  are  also  calculated  at  each  segment. 

The  finite-element  algorithm  is  checked  against  a  test  case  found  in  Pode’s  Tables 
for  Computing  the  Equilibrium  Configuration  of  a  Flexible  Cable  in  a  Uniform  Stream 
[2].  Pode  estimates  the  length  of  a  cable  required  for  a  moored  buoy  to  be  constantly  on 
the  surface  of  the  ocean.  The  buoy  and  cable  are  subject  to  a  constant  current  of  5  kt  in  an 
ocean  depth  of  1097  m.  The  buoy  exhibits  given  vertical  (buoyancy  and  dynamic  lift)  and 
horizontal  (drag)  force  vectors  of  40500  N  and  23100  N,  respectively.  The  cable 
diameter  Z)  =  0.1  m  and  mass  per  unit  length  of  0.4  kg/m  result  in  an  estimated  cable 
specific  gravity  of  4.04  in  13°C  seawater  (seawater  density  pw  =  1026  kg/m  ,  from 
Kinsler  et  al.  [8]).  Pode’s  estimate  of  the  required  cable  length  is  calculated  to  be  2670  m. 


The  MATLAB  algorithm  assumes  the  initial  forces  and  cable  properties  (length, 
specific  gravity)  given  in  Pode’s  example  to  determine  the  shape  of  a  2670  m  cable.  By 
using  the  given  normal  drag  force  on  the  cable  Fd  =  56.9  N/m,  the  cable  drag  coefficient 
Cz)  =  1.5  is  obtained  from  the  rearrangement  of  Munson,  et  al’s  equation  for  Fr,  as  a 
function  of  Cd 


Co  = 


2K 


if) 


whereof  is  the  cross-sectional  area  of  the  cable  of  the  1  m  reference  segment  [9]. 
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Figure  6  shows  the  shape  of  Pode’s  results,  while  Figure  7  graphs  the  MATLAB 
result  for  a  2670  m  long  cable.  For  the  buoy  to  be  at  the  surface  of  a  1097  m  deep  water 
column,  the  MATLAB  algorithm  needs  the  cable  drag  coefficient  to  be  Cd  =  1.6.  The 
cable  shapes  are  similar  and  the  estimated  tilt  angles  relative  to  the  horizontal  at  critical 
points  are  the  same.  For  the  Pode  data,  an  angle  of  cpi  =  171°  relative  to  the  horizontal  is 
found  where  the  buoy  cable  connects  to  the  seafloor  mooring.  Using  the  MATLAB 
algorithm,  this  angle  is  identical  at  171°.  The  angle  where  the  cable  connects  to  the  buoy 
is  the  same  at  (p2  =  120°  in  both  cases.  Note  that  the  scale  of  Figure  7  is  adjusted  to  show 
the  similarity  between  the  figure  found  in  Pode’s  example  and  the  algorithm  results. 


Figure  6.  Drawing  of  Pode’s  example,  5  =  2670  m,  V=  5  kt,  current  flows  right  to  left 

(From  [2]) 
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Figure  7.  Algorithm  results  from  Pode,  5  =  2670  m,  Cd  =  1.6,  Cf  =  0.012,  V=  5  kt, 

current  flows  right  to  left 

B.  SPHERICAL  BODY  FLOW  CHARACTERISTICS 

A  sphere  in  a  horizontal  flow  field  produces  drag  and  buoyancy  forces  that  govern 
the  initial  tensions  on  the  cable.  The  dimensionless  Reynolds  number  (Re)  is  first 
calculated  to  graphically  determine  the  drag  coefficient  Cn,  which  is  proportional  to  the 
horizontal  drag  force  Fq.  The  vertical  force  is  found  by  subtracting  the  buoyancy  force  Fg 
from  the  weight  of  the  sphere  Wspn  and  adding  the  lift  force  Fg.  Re  is  used  to  graphically 
determine  the  lift  coefficient  Ci,  which  is  also  proportional  to  Fl.  The  horizontal  and 
vertical  forces  become  the  horizontal  and  vertical  tension,  th  and  tk  respectively.  Table  2 
lists  the  symbols  and  definitions  used  to  the  find  th  and  ly  imparted  at  the  end  of  the  cable 
by  the  sphere. 
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Table  2.  Spherical  body  symbol  definitions 


Symbol 

Equation  Definitions  and  Units 

g 

Gravity  acceleration  constant  (9.8  m/s  ) 

Re 

Reynolds’s  number 

V 

Flow  velocity  (kt) 

Dsph 

Diameter  of  sphere  (m) 

Pw 

-5 

Seawater  density  (kg/m  ) 

1 

Specific  weight  of  seawater  (N/m  ) 

V 

Kinematic  viscosity  of  seawater  (m  /s) 

p 

Dynamic  viscosity  of  water  (N-s/m) 

Cd 

Drag  coefficient 

Fd 

Drag  force  (N) 

Thi 

Initial  horizontal  cable  tension  (equivalent  to  Fd) 

VsPH 

Volume  of  sphere  =  ttDsph  /6  (m  ) 

PSPH 

Sphere  material  density  (kg/m  ) 

WsPH 

Weight  of  sphere  (N) 

Fb 

Buoyancy  force  (N) 

Cl 

Lift  coefficient 

Fl 

Lift  force  (N) 

TVI 

Initial  vertical  cable  tension  (equivalent  to  Fb  -  Wsph+  Fl) 

Ti 

Resultant  initial  tension  (N) 

0i 

Initial  cable  angle  (deg) 
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To  find  Re,  the  formula  from  Munson  et  al.  [9]  is  used: 

Re  =  V^DIv  (8) 

where  v  =  kinematie  viseosity  of  seawater  is  v  =  ji/pw- 

Assuming  a  temperature  of  13°C,  the  input  values  for  dynamie  viseosity  of  water 

3  2  3 

p  =  1.2*10  N-s/m  [9]  and  seawater  density  pw  =  1026  kg/m  [8]  are  assumed.  The 

6  2 

kinematie  viseosity  for  seawater  is  then  found  to  be  v  =  1.17*10'  m  /s. 

The  drag  eoeffieient  Cd  for  a  sphere  is  found  by  interpolating  the  log-log  graph  in 
Figure  7  of  experimental  data  of  Re  vs.  Cd  from  Constantinescu  and  Squires  [10]. 
Seetions  of  the  graph  are  estimated  by  the  program  using  linear  estimates  of  Cd  from 
segments  of  Re  varying  from  approximately  0  to  4*10^  (within  the  range  of  expeeted  Re 
values)  of  the  Sehliehting  values  on  Figure  8.  More  exaet  solutions  from  other  sourees 
eorrelate  to  Sehliething’s  data  and  are  used  in  the  MATLAB  algorithm.  For  example,  for 
a  Re  between  0  and  1,  Cd  can  be  approximated  from  Stake’s  law  where  Cd  =  24.0/i?e  [9]. 
Also,  the  Cd  plateau  at  5*10^  <  <  2.9*10^  is  estimated  to  be  a  eonstant  value  Cd  = 

0.47  from  examining  experimental  data  in  Hoerner  [4].  A  departure  from  the  trend  near 
“a  eritieal  Reynolds  number  of  RecR  =  3.7*10^”  is  a  steep  dip  that  is  “a  manifestation  of 
the  signifieant  differences  between  laminar  and  turbulent  boundary-layer  separation” 
[10]. 
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Figure  8.  Experimental  data  showing  Co  as  a  function  of  Re  for  uniform  flow  over  a 

sphere  (From  [10]) 

With  the  approximate  Co  known,  the  drag  force  Fd  on  the  sphere  is  calculated  by 
rearranging  Munson,  et  al’s  equation  [9]  for  Cd  as  a  function  of  Fd, 

Fr,  =  (9) 

The  force  in  the  vertical  direction  is  the  buoyancy  force  Fg.  The  buoyancy  force 
of  the  sphere  is  calculated  by  estimating  the  displacement  volume  of  seawater  [9].  The 
specific  weight  of  seawater  at  13°C  is  calculated  as  y  =  1.01*10"^  N/m^  and  is  assumed  to 
be  constant  and  independent  of  depth. 

Fb  =  (10) 

where  y  =  specific  weight  of  seawater  =  p„g  (N/m  ). 
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Another  vertical  force  is  the  lift  generated  by  the  sphere.  Rearranging  Munson,  et 
al.’s  equation  [9],  the  lift  force  Fl  depends  on  the  Cl  and  is  given  as: 


Fl  =  \p.y'^DC,  (11) 

However,  for  a  sphere  at  rest,  the  Cl  for  Re  >  1  is  very  small  as  shown  by 
Buscaglia  [12],  Since  Fl  linearly  depends  on  Cl,  Fl  is  considered  to  be  insignificant  and 
is  neglected  by  the  MATLAB  program. 

With  both  horizontal  and  vertical  forces  known,  the  tension  at  the  free  end  of  the 
cable  system  is  calculated  using  the  horizontal  and  vertical  forces  of  the  sphere.  The 
horizontal  tension  xhi  is  the  same  as  the  drag  force  Fd  (Equation  12).  The  vertical  tension 
xvi  is  the  weight  of  the  sphere  Wsph  (Equation  13)  subtracted  from  the  buoyancy  force  Fg 
(Equation  14).  The  resultant  initial  tension  r,  is  calculated  by  taking  the  magnitude  of  the 
combined  force  vectors  (Equation  15).  The  initial  angle  6i  acting  on  the  object  is  found 
by  taking  the  arctangent  of  the  vertical  and  horizontal  forces  (Equation  16). 


II 

(12) 

^Sph  -  gP^FFVsPH 

(13) 

'^VI  ~  ~  ^SPH 

(14) 

T  ~  HI  ) 

(15) 

e,  =  tan-'(^) 

(16) 
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c. 


STREAMLINE  BODY  FLOW  CHARACTERISTICS 


Similar  to  the  sphere,  the  flow  eharaeteristies  of  a  streamlined  body  produee 
forees  that  govern  the  initial  drag  and  buoyaney  forees  on  the  eable.  The  method  for 
eomputing  the  forees  is  the  same  with  the  differenees  arising  from  the  values  of  Cd  and  A 
for  a  fixed  volume.  The  values  for  Cd  for  a  given  fluid  and  speed  for  the  streamline  body 
are  less  than  a  sphere  of  eomparable  diameter.  As  a  result,  the  streamline  body  has  less 
drag  than  the  sphere  of  similar  volume.  The  dimensional  inputs  for  the  MATLAB 
program  are  the  diameter  of  the  hemispherieal  seetion  (d),  length  of  the  right  eireular 
eylinder  (/),  and  length  of  the  eonie  tail  seetion  (T).  The  overall  length  of  the  streamline 
body  (f)  is  found  by  adding  the  radius  of  the  hemisphere  {d/2),  /,  and  T.  From 
experimental  analysis,  shown  in  Figure  9,  the  drag  eoeffieient  of  the  wetted  surfaee 
{CDwet),  and  thus  Fd,  is  a  funetion  of  the  £/d  ratio  and  Re.  Figure  10  displays  some 
eonfigurations  of  streamline  body  shapes  [4].  For  this  study.  Configuration  B  is 
eonsidered  with  dimensions  given  in  Figure  11. 


Figure  9.  Seleeted  results  on  the  drag  of  rotationally-symmetrie  bodies  (no  eorreetions 

applied)  (From  [4]) 
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BASIC  SHA^e  KIUNO  IN  COM- 
MfRCIAL  AIRRLANC  FUSCUMM 


EXAMPLE  Of’  TXE  SHAPE  OF 
AN  AVERAGE  AIRSHIP  HULL- 


SASIC  SHAPE  Of  TUR0O- 
JFT  ENGINE  NACELLES. 


Figure  10.  Some  examples  of  basie  streamline  body  shapes  in  engineering  applications 

(From  [4]) 


< - >  T 

A 


< - >  I 

< - >  e 

d  s  diameter  of  cylindrical  part  of  streamline  body  (m) 

€  =  length  of  streamline  body  (m) 

/  =  length  of  right  circular  cylinder  section  of  streamline  body  (m) 

T  s  length  of  conic  tail  section  of  streamline  body  (m) 

Figure  1 1 .  Dimensions  used  for  streamline  body 
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Table  3.  Streamline  body  symbol  definitions 


Symbol 

Equation  Definitions  and  Units 

g 

Gravity  aeeeleration  eonstant  (9.8  m/s  ) 

Re 

Reynolds  number 

V 

Flow  veloeity  (kt) 

d 

Diameter  of  eylindrieal  part  of  streamlined  body  (m) 

€ 

Length  of  streamline  body  (m) 

1 

Length  of  the  right  eireular  eylinder  seetion  of  streamline  body  (m) 

T 

Length  of  eonie  tail  eone  seetion  of  streamline  body  (m) 

Cjlam 

Laminar  flow  skin  frietion  drag  eoeffieient 

Cf 

Turbulent  flow  skin  frietion  drag  eoeffieient 

Cowet 

Drag  eoeffieient  for  streamline  body 

Fd 

Drag  foree  (N) 

t^hi 

Initial  horizontal  eable  tension  (equivalent  to  Fq) 

VsB 

Volume  of  streamline  body  (m  ) 

PSB 

Streamline  body  material  density  (kg/m  ) 

WsB 

Weight  of  streamline  body  (N) 

Fb 

Buoyaney  foree  (N) 

Cl 

Lift  eoeffieient 

Fl 

Lift  foree  (N) 

TVI 

Initial  vertieal  eable  tension  (equivalent  to  Fb  -  Wsph+  Fl) 

Ti 

Resultant  initial  tension  (N) 

Oi 

Initial  eable  angle  (deg) 

20 


The  method  for  finding  the  drag  forces  and  buoyancy  forces  is  the  same  as  that  of 
the  sphere.  The  Re  of  the  streamline  body  is  obtained  from  Equation  8  except  that  d,  the 
diameter  of  the  cylindrical  section,  is  substituted  for  D.  With  Re  known,  the  drag 
coefficient  for  the  streamline  body  Cowet  is  found.  Figure  13  displays  the  values  for  Cowet 
relating  to  Re  for  a  streamline  body  [4].  The  program  models  Cowet  as  a  function  of  Re  in 
two  conditions:  Re  <  10^  and  Re  >  10^.  It  is  assumed  that  the  streamline  body  shown  in 
Figure  12  follows  the  same  as  those  of  Configuration  B  with  the  same  length/diameter 
ratios. 

For  Re  of  lO"^  <  Re  <  10^,  the  flow  is  laminar.  The  Blasius  solution  for  Cjiam 
(laminar  skin-friction  coefficient)  for  flow  over  a  fiat  plate  characterizes  the  drag.  To 
account  for  the  size  of  the  streamline  body,  Cowet  is  a  function  of  Cjiam  and  the  diameter- 
to-length  ratio  {d/€).  Note  that  for  an  extremely  large  d/i,  the  streamline  body  resembles  a 
fiat  plate  and  Crhvet  =  Cjiam-  For  i?e<10"^,  the  MATFAB  program  assumes  that  the  Cowet 
follows  the  Blasius  solution  for  laminar  flow. 


Blasius  solution  for  flow  over  a  fiat  plate  for  low  Re  [9]  gives 


C 


flam 


1.328 

fRe 


(17) 


The  drag  coefficient  for  streamline  body  with  laminar  flow  [4]  is 


+  (18) 

At  higher  (10^  <  <  10'°),  the  flow  is  turbulent.  Hoemer  uses  data  from  the 

experimental  data  from  the  Schoenherr  line  for  flow  over  a  surface  (Figure  12)  to 
compute  Equation  19  for  the  skin-friction  drag  coefficient  C/.  To  account  for  the  effect  of 
streamline  body  dimensions,  Hoerner  developed  Equation  20  to  find  Cowet  using 
experimental  data  in  Figure  13  turbulent  flow  over  a  streamline  body. 
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Figure  12.  Average  or  total  skin-frietion  drag  eoeffieient  of  smooth  and  plane  surfaees 

(in  ineompressible  flow)  in  air  and  in  water  (From  [4]) 


The  turbulent  skin-friction  drag  coefficient  using  the  Schoenherr  line  is 
approximated  by 


'  [3.461og(i?e)-5.6]" 


(19) 


•  rUBRKAVB  3H1PK  BO.  1  WITH  dA  -  ^  6  %  (21  .g) 

«  AHC-iroB  "B*BODT  WITH  dA  “  20  (21  .h) 

X  BACA-VDT  WITH  dA»  12*0  16  <  (21 , a) 

d  BAOA-VDT  WITH  dA  -  17  TO  22*  (21  ,l)AHIl(22,k) 

•  BACA  AIRSHIP  HULL  WITH  dA  -  10  *  (21  ,d) 

•  BSRLIB  IB  T0WIH9  lAHK  WITH  d/1  -  (22, o) 

•  HARBURO  IB  TOWIHO  TABK  WITH  dA  -  17  it  (22.1>) 


Figure  13.  Drag  of  streamline  bodies,  tested  in  turbulent  wind  tunnels  or  with  turbulence 

simulation  (From  [4]) 
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As  with  the  sphere,  Cowet  is  found  and  inputted  into  Equation  9  to  ealeulate  Fd. 
The  buoyant  foree  is  found  from  Equation  10  exeept  now  the  volume  is  that  of  the 
streamline  body  using  Equation  21. 


Fsb  =  volume  of  hemisphere  +  volume  of  right  eireular  eylinder  +  volume  of  eone 

Vss  =  -d^  +  -d^L  +  —d^T 
6  4  12 


(21) 


Eor  simplieity,  the  length  of  the  tail  eone  T  is  set  to  be  the  same  as  the  diameter  of 
the  eylinder  d.  Fi  is  assumed  to  be  negligible  sinee  the  streamline  body  is  assumed  to 
have  uniform  pressure  (hydrostatie)  at  all  points  around  the  body  and  a  zero  angle  of 
attaek  (no  lift  is  generated).  The  tension  r,  and  angle  9i  on  the  modem  produeed  by  the 
streamline  body  is  ealeulated  with  Equations  15  and  16. 
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IV.  CASE  STUDY-THE  MOORED  CABLE 


The  first  case  study  is  a  modem  moored  in  the  deep  ocean.  The  ocean  depth  is 
assumed  to  be  3000  m.  The  desired  depth  of  the  modem  is  1000  m  with  the  cable  being 
approximately  2000  m  long  (see  Figure  14). 

Seasurface 


Syntactic 
Foam  Float 


Cylindrical 


Figure  14.  Moored  case  study  with  spherical  float 

BZ-26  syntactic  foam  is  used  as  flotation  material.  The  cable  is  chiefly  comprised 
of  Nitronic  50B  Stainless  Steel  and  has  a  net  diameter  of  16  mm  (5/8  in).  The  cable  is 
typically  made  of  many  strands  of  material  braided  into  a  cylindrical  shape,  and  the  net 
diameter  is  the  diameter  of  all  the  strands  combined.  It  is  expected  that  the  net  Yield 
Strength  (<775)  of  the  cable  is  greater  than  that  of  a  solid  material’s  ars  with  the  individual 
strands  acting  together.  However,  the  ays  assumed  in  these  case  studies  is  that  of  the  solid 
material  itself  to  represent  a  lower  bound  for  material  failure.  It  is  assumed  the  cable  is 
fabricated  in  such  a  way  that  it  is  flexible  and  supports  a  wide  range  of  bend  angles.  At 
this  diameter,  with  a  ays  of  380  MPa  for  Nitronic  50B,  the  maximum  allowable  tension  in 

the  cable  (Zmax)  is  75.2  kN  before  plastic  deformation. 
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^max  =  (22) 

The  determination  of  the  appropriate  step  size  ds  for  the  finite-element 
representation  is  important  to  accurately  estimate  array  position.  As  the  step  size 
decreases,  the  algorithm  conducts  more  calculations  and  requires  more  computation  time. 
There  is  a  point  of  diminishing  returns  at  which  the  step  size  is  small  enough  that  using  a 
smaller  step  size  does  not  appreciably  improve  the  results.  The  minimum  step  size  that 
can  be  used  with  this  MATLAB  program  is  =  1  m.  To  find  the  most  efficient  step  size, 
the  case  of  a  2.4  m  diameter  spherical  float  in  a  current  of  3  kt  is  presented.  Using  BZ-26 
syntactic  foam,  the  buoyancy  force  is  72.8  kN.  The  catenary  is  plotted  for  step  sizes  of 
100,  10,  and  1  m.  Figure  15  shows  that  at  a  step  size  of  =  1  m  the  catenary  has 
converged.  Therefore,  for  the  remaining  analysis  for  the  2000  m  cable  uses  step  size 
ds=  \  m. 
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Figure  15.  Step  size  variations  for  2.4  m  diameter  spherical  float,  Nitronic  50B  cable,  D 
=  16  mm,  s  =  2000  m,  F  =  3  kt,  current  flows  right  to  left 
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The  volume  of  the  flotation  material  relative  to  the  eable  properties  impaets  the 
vertieal  displacement.  For  a  given  cable  length,  if  the  buoyancy  force  is  not  large  enough, 
the  weight  of  the  cable  can  sink  the  modem  below  the  desired  operating  depth.  The 
weight  of  a  2000  m,  16  mm  diameter  Nitronic  SOB  cable  is  3090,  thus  requiring  a  large 
buoyancy  force  to  suspend  the  cable.  The  algorithm  indicates  that  the  vertical  forces  at 
the  end  of  the  cable  are  inadequate  when  the  cable  crosses  the  imaginary  seafloor 
boundary.  This  is  shown  in  Figure  16  when  the  cables  associated  with  spherical  diameters 
of  2.1  and  2.0  m  intersect  the  seafloor.  Although  such  interaction  with  the  ground  would 
change  the  forces  acting  on  the  cable,  the  computations  in  this  algorithm  do  not  account 
for  the  seafloor. 


2>0  PoMion  of  Moo(>d  Cab)* 


Figure  16.  Spherical  float  diameter  variations,  Nitronic  SOB  cable,  D=  \6  mm,  s  =  2000 

m,  V  =  3  kt,  current  flows  right  to  left 
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The  shape  of  the  float  strongly  influences  the  amount  of  horizontal  drag  Fd 
imposed  on  the  end  of  the  array.  The  sphere  has  a  higher  Cd  and  a  higher  cross-sectional 
area  than  the  streamline  body  of  the  same  volume.  The  volume  is  chosen  to  elevate  the 
modem  at  a  prescribed  depth  across  a  given  current  range.  For  this  case  study,  to  get  the 
same  buoyancy  force  as  a  2.3  m  diameter  spherical  float,  a  streamline  body  with  1.07  m 
diameter  and  6.96  m  length  is  specified.  With  a  current  of  3  kt,  the  sphere  Cd  is  0.232  and 
the  streamline  body  Cd  is  0.004.  The  Fd  on  the  sphere  is  1.28  kN  and  on  the  streamline 
body  is  9.6  N.  Figure  17  compares  the  cable  shapes  produced  by  the  sphere  and 
streamline  body.  Note  that  with  less  drag  force  on  the  streamline  body,  the  horizontal 
position  of  the  array  changes  less  over  the  depth  range  and  also  the  maximum  tension  on 
the  cable  is  less  with  x sphere  =  45.2  kN  and  xsb  =  45.1  kN.  Both  tensions  are  less  than  Xmax 
by  nearly  a  factor  of  1.6.  Note  that  the  shapes  of  the  cable  for  each  shape  are  similar.  For 
such  a  long  length,  the  drag  and  lift  forces  acting  on  the  cable  dominate  the  cable 
displacement. 
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Figure  17.  Comparison  of  sphere  vs.  streamline  body  float,  Nitronic  50B  cable,  D=  \6 

mm,  s  =  2000  m,  V  =  3  kt,  current  flows  right  to  left 
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For  further  analysis,  the  streamline  body  is  ehosen.  As  the  eurrent  is  inereased 
from  0  to  3  kt,  the  shape  of  the  eable  eatenary  ehanges  dramatieally.  At  low  speeds,  the 
eable  is  nearly  vertical.  As  current  speed  increases,  Fd  increases  to  cause  more 
horizontal  displacement  resulting  in  a  modem  depth  deeper  than  the  desired  1000  m. 
Figure  18  shows  cable  shape  for  the  array  for  a  streamline  body  with  diameter  d  =  1.07 
m  diameter  and  length  i  =  6.96  m. 


2-0  PoMixi  of  Mootod 
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Figure  18.  Comparison  of  cable  shapes  with  streamline  body  float,  Nitronic  50B  cable,  V 
varies  0-3  kt,  Z)  =  16  mm,  s  =  2000  m,  current  flows  right  to  left 
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The  cable  material  is  varied  when  a  “high-strength,  light  weight  polyethylene 
fiber”  is  substituted  for  Nitronic  SOB  [12],  Honeywell’s  Spectra  fiber  1000,  whose 
properties  are  in  Appendix  C,  has  a  high  resistance  to  chemicals,  water,  ultraviolet  light 
and  flexural  fatigue.  The  product  has  good  resistance  to  abrasion  and  is  used  in  marine 
applications  such  as  marine  lines  as  well  as  fishing  lines  and  commercial  fishing  nets. 

-5 

The  product  density  is  970  kg/m  ,  resulting  in  a  specific  gravity  of  0.95  in  13°C  seawater. 
For  comparison  of  the  effects  of  a  lighter  material  on  cable  shape,  a  2000-m  15.88-mm 
diameter  Spectra  1000  (1300  denier)  line  is  modeled  against  a  16-mm  diameter  Nitronic 
SOB  cable  at  3  kt.  The  maximum  tension  found  in  the  Spectra  1000  line  is  45.4  kN  while 
Tmax  (based  on  the  products’  Ultimate  Tensile  Strength  of  3.00  GPa)  is  594kN.  As 
expected,  the  buoyant  Spectra  has  less  downward  vertical  displacement  than  the  heavier 
Nitronic  5 OB. 
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Figure  19.  Comparison  of  Nitronic  50B  vs.  Spectra  1000  cables  with  streamline  body 
float,  D=  \6  mm,  s  =  2000  m,  V  =  3  kt,  current  flows  right  to  left 
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The  diameter  of  the  eable  ean  also  be  adjusted  to  find  the  minimum  eable 
diameter  before  eable  failure.  The  eable  diameter  in  Figure  20  were  varied  to  find  when 
breakage  would  oeeur.  By  only  ehanging  the  cable  diameter  D,  it  was  found  that  the 
cable  size  could  be  reduced  to  12.7  mm  (1/2  in).  The  maximum  tension  found  in  the  cable 
of  45.3  kN  nearly  exceeds  the  calculated  Xmax  of  48.1  kN. 

The  smaller  cable  diameter  produces  less  cable  drag  thus  less  horizontal  force, 
resulting  in  a  more  vertical  cable.  The  tradeoff  with  reduced  cable  diameter  is  that  there  is 
less  material  and  inherently  less  allowable  strength.  Since  this  is  a  steady-state  analysis, 
snap  loads  should  be  considered  when  actually  specifying  the  cable  dimensions. 
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Figure  20.  Comparison  of  Spectra  1000  cable  with  varying  cable  diameter  D,  streamline 
body  float,  5  =  2000  m,  V  =  3  kt,  current  flows  right  to  left 
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V.  CASE  STUDY-THE  TOWED  CABLE 


The  second  case  study  is  an  acoustic  modem  suspended  from  a  surface  buoy  in 
the  ocean.  The  desired  depth  of  the  modem  is  30  m  with  the  cable  scope  being 
approximately  30  m.  The  cable  is  chiefly  comprised  of  Nitronic  50B  stainless  steel  with 
an  enclosed  electrical  connection  between  the  buoy  and  the  modem.  The  cable  has  a  net 
diameter  of  6.4  mm  (1/4  in).  At  this  diameter,  the  maximum  allowable  tension  in  the 
cable  is  found  to  be  imax  =  12.0  kN  (using  Equation  22).  The  same  assumptions  as  with 
the  cable  in  the  moored  case  study  are  made  here. 
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Figure  21 .  Towed  case  study  with  spherical  ballast 


The  same  analysis  done  with  the  moored  cable  can  be  performed  on  the  towed 
cable.  First,  the  most  appropriate  step  size  ds  for  the  finite-element  algorithm  is  found  by 
parametric  analysis  shown  in  Figure  22.  Convergence  of  cable  shape  occurs  with  step 
size  ds  =  0.\  m.  For  the  remainder  of  this  case  study  analysis,  a  step  size  of  =  0.1  m  is 
used. 
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Figure  22.  Step  size  variations  for  a  0.2  m  diameter  spherical  ballast,  Nitronic  50B  cable, 
D  =  6.4  mm,  5  =  30  m,  F  =  3  kt,  current  flows  right  to  left 

The  size  of  the  ballast  material  relative  to  the  cable  impacts  the  vertical 
displacement.  A  larger  ballast  of  similar  material  pulls  the  pulls  the  end  of  the  cable  down 
more.  If  the  ballast  is  too  small,  it  may  not  be  enough  to  maintain  the  modem  at  the 
desired  depth.  The  algorithm  indicates  that  the  vertical  forces  at  the  end  of  the  cable  are 
inadequate  when  the  cable  crosses  the  imaginary  seasurface  boundary.  When  the  cable  is 
light,  and  the  size  of  the  ballast  is  very  small,  the  cable  will  appear  to  be  above  the  sea 
surface.  Although  the  cable  cannot  actually  breach  the  upper  boundary  of  the  sea  surface, 
the  MATLAB  algorithm  does  not  recognize  a  boundary  violation. 
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As  previously  determined,  the  sphere  has  a  higher  Cd  and  A  than  the  streamlined 
body.  A  0.2  m  diameter  spherical  ballast  has  the  same  weighting  characteristics  as  a 
streamline  body  of  0.09  m  diameter  and  0.58  m  length.  For  a  current  of  3  kt,  the  drag 
coefficient  of  the  sphere  is  Cd  =  0.47  and  drag  force  is  Fd  =  18.0  N.  The  streamline  body 
values  are  Cd  =  0.008  and  Fd  =  0.12  N.  Figure  23  compares  the  cable  shape  arising  from 
a  sphere  and  streamline  body.  Again,  with  less  drag  on  the  streamline  body,  the 
horizontal  displacement  of  the  cable  is  less.  The  maximum  tension  on  the  cable  is  nearly 
the  same  with  xsphere  =  341  N  and  tsb  =  342  N.  These  tensions  are  less  than  Xmax- 
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Figure  23.  Comparison  of  sphere  vs.  streamline  body  ballast,  Nitronic  SOB  cable,  D  =  6.4 

mm,  5  =  30  m,  F  =  3  kt,  current  flows  right  to  left 


35 


The  streamline  body  ballasted  cable  shape  is  predicted  for  currents  of  0  to  3  kt. 
The  cable  behavior  is  similar  to  the  moored  cable  where  at  lower  speeds,  the  cable  is 
nearly  vertical  and  at  higher  speeds,  depth  decreases  while  horizontal  displacement 
increases.  Figure  24  shows  the  shape  of  the  cable  for  a  streamline  body  with  d  =  0.09  m 
and  €  =  0.58  m. 
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Figure  24.  Comparison  of  cable  shapes  of  streamline  body  ballast  as  V  varies  0-3  kt,  D  = 

6.4  mm,  5  =  30  m,  current  flows  right  to  left 
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In  Figure  25,  Spectra  1000  is  substituted  for  Nitronic  50B  using  the  streamline 
body.  With  the  lighter  Spectra  1000,  the  cable  is  less  vertical  as  the  buoyancy  forces  act 
in  an  upward  direction.  For  the  given  cable  diameter  D  =  6.4  mm,  the  maximum 
allowable  tension  in  the  cable  with  the  Spectra  1000  is  found  to  be  95.0  kN,  compared  to 
12.0  kN  with  the  Nitronic  50B.  The  analysis  for  the  minimum  diameter  cable  needed 
under  these  conditions  can  also  be  done. 


2-D  P»MionofTo«fMjCat4i 


Figure  25.  Comparison  of  Nitronic  50B  vs.  Spectra  1000  cables  with  streamline  body 
ballast,  D  =  6.4  mm,  5  =  30  m,  F  =  3  kt,  current  flows  right  to  left 
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VI.  CONCLUSIONS 


A  horizontal  current  acting  on  a  vertical  cable  imparts  a  drag  force  that  displaces 
the  cable.  To  find  the  overall  cable  shape,  first  the  tension  and  initial  angle  imparted  by 
the  vertical  and  horizontal  forces  acting  on  the  ballast/fioat  shape  are  found.  With  the 
initial  conditions  known,  the  drag  and  lift  acting  on  the  cylindrical  modem  and  uniform 
cable  are  found  and  accumulated  over  the  length  of  the  array  to  the  fixed  point.  With  the 
array  tensions  and  angles  known,  the  cable  shape  is  determined. 

Depending  on  the  length  of  the  array,  the  ballast/float  shape  may  or  may  not 
significantly  influence  the  array  shape.  A  streamline  body  compared  to  a  sphere  of  equal 
volume  has  less  initial  drag  force.  However,  when  the  cable  is  very  long  (as  in  the  2000- 
m  case)  the  position  of  the  array  is  more  dependent  on  drag  of  the  cable.  At  a  shorter 
cable  scope,  the  forces  on  the  fioat/ballast  body  are  more  dominant.  In  either  case,  it  is 
recommended  that  a  streamline  body  shape  be  used  as  the  initial  drag  forces  are  less. 

Design  of  a  tethered  modem  requires  knowing  the  desired  operational 
requirements  of  the  system.  The  system  tolerances  dictate  the  cable  and  ballast/fioat 
material  type  and  size  and  the  designer  should  understand  the  trade-offs.  For  example, 
while  a  larger  float  shape  may  add  buoyancy,  the  additional  drag  may  actually  move  it 
farther  from  the  desired  position.  Using  this  MATLAB  algorithm,  the  designer  can 
analyze  these  factors  and  specify  a  design  to  meet  the  operational  requirements. 
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APPENDIX  A.  MATLABCODE 


%  CABLE  ESTIMATOR  is  a  program  to  compute  the  initial  value  problem  for  a 
%  first  order  system  of  differential  equations  describing  the  shape  of  an 
%  array  system.  The  array  can  either  be  towed  or  moored  on  the  bottom  of 
%  the  ocean  and  has  a  eylindrical  acoustic  modem  attached  at  the  end.  The 
%  program  allows  the  user  to  consider  two  types  of  ballast  or  float  shapes 
%  attached  at  the  end  of  the  modem  -  a  sphere  or  streamline  body.  The  user 
%  can  also  input  in  arbitrary  initial  tensions  for  towed  or  moored  cases. 

%  The  program  allows  the  user  to  change  the  material  properties  and 
%  dimensions  of  the  cable.  The  outputs  of  the  program  are  tension,  the 
%  angle  relative  to  the  horizontal,  and  x,  y  position  on  a  two-dimensional  plot. 

%  Symbol  Definitions 

%  A  =  eross-sectional  area  of  eable  at  a  point  (m^2) 

%  DENC  =  eable  density  (kg/m^3) 

%  DENW  =  water  density  (kg/m^3) 

%  thetadeg  =  the  angle  between  the  horizontal  and  local  tangent  at  a  point 
%  on  the  cable  (deg)  (positive  x  axis  is  0  deg) 

%  V  =  towing  velocity  (or  current)  (kt) 

%  DC  =  diameter  of  cable  (m) 

%  tauN  =  cable  tension  (N) 

%  s  =  cable  length  (m;  s  =  0  is  end  of  array) 

%  Co  =  basic  drag  coefficient 

%  Cf  =  drag/lift  eoefficient  eorreetion  factor  (Cf  =  pi  *  Cp;  Cp  is  the 
%  classical  drag  coefficient;  See  Hoerner  pg.  3-11  "Cross-Elow  Prineiple") 

%  Cl  =  correeted  lift  coefficient 

%  Cd  =  corrected  drag  eoefficient 

%  g  =  gravity  acceleration  constant  (m/s^2) 

%  1  m  =  3.2808  ft  =  39.37  in 
%  1  kg  =  2.2046  lb 
%  1  Ibf  =  4.45  N 

clear  all 

%Program  Inputs 
V  =  3.0;  %knots 

Shape  =  2;  %  Shape  =  1  =>  Sphere,  Shape  =  2  =>  Streamline  body, 

%  Shape  =  3  =>  No  attached  body  or  sensor  (towed  cable), 

%  Shape  =  4  =>  No  attaehed  body  or  sensor  (moored  eable) 

%  If  Shape  =  3  or  4,  input  THI  and  TVI  on  lines  283/284  for  arbitrary 
%  initial  horizontal  and  vertical  forces. 
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DSph  =  2.3;  %  Diameter  of  sphere  (m) 

%  Float  2.3  m 
%  Ballast  0.2  m 

DSB  =  1.071091;  %  Maximum  diameter  of  streamline  body  (m) 

%  Float  1.071091  m  (Vol  2.3  m  sphere) 

%  Ballast  0.08933  m  (Vol  0.2  m  sphere) 

LSB  =  6.962093;  %  Length  of  cylindrical  part  of  streamline  body  (m) 

%  Float  6.962093  m  (Vol  2.3  m  sphere) 

%  Ballast  0.58069  m  (Vol  0.2  m  sphere) 

TSB  =  1.071091;  %  Length  of  conical  section  (tail)  of  streamline  body  (m) 

%  Float  1.071091  m  (Vol  2.3  m  sphere) 

%  Ballast  0.08933  m  (Vol  0.2  m  sphere) 

Density  =  1;  %  Density  of  the  shape,  DEN  =  1  =>  DENF  (float), 

%  DEN  =  2  =>  DENB  (ballast) 

E  =  2000;  %  Cable  length  (m)  Pode  example  2760  m 

DC  =  0.015875;  %  Cable  diameter  (m)  Moored  0.015875  Towed  0.00635  Pode 
0.0111125 

ds  =  1;  %  Step  interval  (m)  limit  is  0.1 

SGC  =  0.95;  %  Specific  gravity  of  cable  in  sw 
%  Nitronic  50B  7.6096  Pode  4.04  SpectralOOO  0.95 

%  Common  Terms 
g  =  9.80665;  %  m/s'^2; 

V  =  0.51444  *  V;  %  Conversion  from  knots  to  m/s 

DENB  =  7807.4;  %  Density  of  Nitronic  50B  (kg/m^3)  for  ballast  material  (Calculated) 
DENE  =  390;  %  Density  of  BZ-26  foam  (kg/m^3)  for  float  material  (Spec  Sheet) 
DENW  =  1026;  %  Density  of  seawater  @  13C  kg/m^3  (KECS) 

Dyna  =  1.17  *  10^-3;  %  Dynamic  viscosity  for  water,  "u"  (mu),  Ns/m^2, 

%  water  @  13C,  Munson  p.759) 

Kine  =  Dyna  /  DENW;  %  Kinematic  viscosity  for  seawater  ,  "v"  (nu),  m^2/s, 

%  Munson  p.  19)  (for  sw  L1404e-006) 

gammasw  =  DENW  *  g;  %  Specific  weight  of  seawater  at  13C,  N/m^3,  assume 
%  constant  and  independent  of  depth  (see  Munson  example,  p.71) 

YS  =  380*10^6;  %  Yield  strength 
%  Nitronic  50B  stainless  steel  380  MPa 

%  Spectra  1000  3000  MPa  (Use  ultimate  tensile  strength  as  approximation) 

MaxCTau  =  pi  *  (DC/2)^2  *  YS; 
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if  Density  ==  1 
DEN  =  DENE; 

else  if  Density  ==  2 
DEN  =  DENB; 

end 

end 

%  Elniform  cable  parameters 

DENC  =  SGC  *  DENW;  %  Density  of  cable  (kg/m'^3) 

Co  =  1.1;  %  Value  of  Co  =  1.1  from  Carlson  and  Smith  (p.3) 

%  Erom  experimental  measurements.  To  match  Pode  example,  Co  is  about  1.6. 
Cf  =  0.012;  %  Value  of  Co  =  0.012  from  Carlson  and  Smith  (p.3) 
n  =  E/ds;  %  Number  of  steps  to  calculate  cable  parameters 

%  Cylindrical  sensor  modem  and  weight/float  are  used  to  calculate  the 
%  initial  tension  at  the  end  of  the  cable 

%  Cylindrical  acoustic  modem  dimensions.  These  values  remain  constant. 

DS  =  0.1016;  %  Nominal  acosutic  modem  diameter  (m)  (4  in) 

LS  =  0.2794;  %  Nominal  acoustic  modem  length  (m)  (11  in) 

VOES  =  pi  *  DS^2  *  ES  /  4;  %  Volume  (m^3) 

MS  =  3;  %  Mass  of  acoustic  modem  in  air  (kg) 

DENS  =  MS  /  VOES;  %  Acoustic  modem  density  (kg/m^3) 

AS  =  pi  *  DS  *  DS  /  4;  %  Cross-sectional  area  of  acoustic  modem 
CoS  =  1.1;  %  Same  value  as  cable 
CIB  =  0.012;  %  Same  value  as  cable 

%  Initial  veritcal  and  horizontal  tension  acting  on  the  acosutic  modem  are 
%  found  by  calculating  the  drag  force  (horizontal)  and  buoyancy  force 
%  (vertical)caused  by  the  ballast/float. 

%  Assume  that  the  ballast/float  shapes  cause  negligible  lift. 

%  Assume  the  distance  between  the  ballast/float  and  modem  is  negligible. 

if  Shape  ==  1  %  Sphere 

ReSph  =  V  *  DSph  /  Kine;  %  Renoylds  number  for  sphere 

%  Approximations  of  Cd  for  a  sphere  given  the  expected  Renyolds 
%  number  with  problem  parameters  from  Constantinescu  and  Squires  p.2 
%  (Schlichting)  and  Hoerner  p.3-8. 

if  ReSph>=0  &&  ReSph<l 
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CdSph  =  24.0/ReSph; 
end 

%  Low  Re  approximation  (Munson  p.524)  (Stoke's  Law) 

ifReSph>=l  &&  ReSph<4 
CdSph  =  -5.3333  *  ReSph  +  29.3333; 
end 

if  ReSph>=4  &&  ReSph<10^1 
CdSph  =  -0.5833  *  ReSph  +  10.333; 
end 

if  ReSph>=10^1  &&  ReSph<5*10^1 
CdSph  =  -0.075  *  ReSph  +  5.25; 
end 

if  ReSph>=5*10^1  &&  ReSph<3*10^2 
CdSph  =  - 1.6*1 0^-3  *  ReSph +  0.52; 
end 

if  ReSph>=3*10^2  &&  ReSph<2*10^3 
CdSph  =  -1.4*  10^-4  *  ReSph +  0.745; 
end 

if  ReSph>=2*10^3  &&  ReSph<2.9*10^5 
CdSph  =  0.47; 
end 

%  Plateau  estimation  (Hoerner  p.3-8) 

if  ReSph>=2.9*10^5  &&  ReSph<3.7*10^5 
CdSph  = -4.75*10^-6  *  ReSph  +  1.8475; 
end 

if  ReSph>=3.7*10^5  &&  ReSph<7*10^5 
CdSph  =  2.3256*10^-8  *  ReSph  +  0.0814; 
end 

if  ReSph>=7*10^5  &&  ReSph<2*10^6 
CdSph  =  7.6923*10^-8  *  ReSph  +  0.04615; 
end 

if  ReSph>=2*10^6  &&  ReSph<4*10^6 
CdSph  =  2.5*  10^-8  *  ReSph +  0.15; 
end 
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ifReSph>4*10^6 
CdSph  =  0.25; 
end 

DragSph  =  0.5  *  DENW  *  V^2  *  pi/4  *  DSph^2  *  CdSph;  %  Drag  force  caused  by 
sphere 

VSph  =  pi  *  DSph^S  /  6;  %  Volume  of  sphere 

FBSph  =  gammasw  *  VSph;  %  N  (See  Munson  book  example  p.71) 

WSph  =  g*  DEN  *  VSph;  %  N 

TVI  =  FBSph  -  WSph;  %  N 
THI  =  -DragSph;  %  N 

%  Conditions  imposed  by  ballast  /  float  on  the  acoustic  modem. 
tauB  =  sqrt(TVD2  +  THD2); 
thetaB  =  atan(TVI/THI); 

%  Calculate  the  lift  and  drag  coefficients  on  the  acoustic  modem. 

%  "sign"  is  the  Signum  function  for  MATLAB  (theta  >  0  =>  1,  theta  =  0  =>  0, 

%  theta  <  0  =>  -1 


CIS  =  CoS  *  (sin(thetaB))^2  *  cos(thetaB)  *  sign( thetaB); 

CdS  =  CoS  *  (sin(thetaB))^3  *  sign(thetaB)  +  CfS; 

dfJauS  =  ((g  *  AS  *  (DENS  -  DENW)  *  sin(thetaB))  +  (1/2  *  DENW  *  (V^2)  *  DS 
*  (CdS  *  cos(thetaB)  -  CIS  *  sin(thetaB))))  *  ES; 

tauS  =  tauB  +  dftauS; 

dfjhetaS  =  ((1  /  tauS)  *  (g  *  AS  *  (DENS  -  DENW)  *  cos(thetaB)  -  1/2  *  DENW  * 
(VA2)  DS  *  (CdS  *  sin(thetaB)  +  CIS  *  cos(thetaB))))  *  ES; 

thetaS  =  thetaB  +  df_thetaS; 

TH(1)  =  tauS  *  cos(thetaS); 

TV(1)  =  tauS  *  sin(thetaS); 

x(l)  =  0; 
y(l)  =  0; 

%  Both  x(l)  and  y(l)  provide  normalization  for  the  array.  The  starting 
%  point  of  (0,0)  is  where  the  modem  connects  to  the  cable. 

%  Assume  the  dimensions  of  the  ballast  and  modem  are  negligible  for 
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%  cable  x,y  displacement  calculations. 


theta(l)  =  thetaS; 
tau(l)  =  tauS; 

else  if  Shape  ==  2  %  Streamline  body 
ReSB  =  V  *  DSB  /  Kine; 

LOSB  =  DSB/2  +  LSB  +  TSB;  %  Overall  length  of  streamline  body 
if  ReSB<10^5 

Cflam  =  1 .328  /  sqrt(ReSB);  %  Laminar  skin  friction  drag  coeffcient  (Hoemer 

p.2-4) 

Cdwet  =  Cflam  *  (1  +  ((DSB/LOSB^l.S))  +  (0.1  l*((DSB/LOSBr2)); 

%  This  equation  is  for  subcritical  Re  (lOM  to  10^5).  (Hoerner  p.6-16  Eqn  24) 

%  Assumption  that  this  is  viable  for  Re  <  10^5. 
end 

if  ReSB>=10^5  &&  ReSB<=10^9 

Cf  =  (3.46*logl0(ReSB)-5.6)^-2;  %  (Hoerner  p. 2-5.  Approximation  of 
Schoenherr 

%  line  on  p.2-6,  Figure  5) 

Cdwet  =  Cf  *  (1  +  1.5*(DSB/LOSB)^1.5  +  7*(DSB/LOSB)^3);  %This  equation 
is 

%  for  Re>10^5.  (Hoemer  p.6-17) 
end 

DragSB  =  DENW  *  V^2  *  pi/4  *  DSB^2  *  Cdwet;  %  Drag  force  caused  by  streamline 
body 

VSB  =  (pi/6)  *  DSB^3  +  (pi/4)  *  DSB^2  *  LSB  +  (pi/12)  *  DSB^2  *  TSB;  %  (m^3) 

%  Approximates  volume  of  streamline  body  as  hemisphere  +  cylinder  +  cone. 

FBSB  =  gammasw  *  VSB;  %  N 
WSB  =  g  *  DEN  *  VSB;  %  N 

TVI  =  FBSB  -  WSB;  %  N 
THI  =  -DragSB;  %  N 

%  Conditions  imposed  by  ballast  /  float  on  the  acoustic  modem. 
tauB  =  sqrt(TVH2  +  THH2); 
thetaB  =  atan(TVI/THI); 

%  Calculate  the  lift  and  drag  coefficients  on  the  acoustic  modem. 

%  "sign"  is  the  Signum  function  for  MATLAB  (theta  >  0  =>  1,  theta  =  0  =>  0, 

%  theta  <  0  =>  -1 

CIS  =  CoS  *  (sin(thetaB))^2  *  cos(thetaB)  *  sign( thetaB); 

CdS  =  CoS  *  (sin(thetaB))^3  *  sign(thetaB)  +  CfS; 
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dfJauS  =  ((g  *  AS  *  (DENS  -  DENW)  *  sin(thetaB))  +  (1/2  *  DENW  *  (V^2)  *  DS 
*  (CdS  *  cos(thetaB)  -  CIS  *  sin(thetaB))))  *  ES; 
tauS  =  tauB  +  dftauS; 

dfjhetaS  =  ((1  /  tauS)  *  (g  *  AS  *  (DENS  -  DENW)  *  cos(thetaB)  -  1/2  *  DENW  * 
(VA2)  DS  *  (CdS  *  sin(thetaB)  +  CIS  *  cos(thetaB))))  *  ES; 
thetaS  =  thetaB  +  df_thetaS; 

TH(1)  =  tauS  *  cos(thetaS); 

TV(1)  =  tauS  *  sin(thetaS); 

x(l)  =  0; 
y(l)  =  0; 

%  Both  x(l)  and  y(l)  provide  normalization  for  the  array.  The  starting 
%  point  of  (0,0)  is  where  the  modem  conneets  to  the  cable. 

%  Assume  the  dimensions  of  the  ballast  and  modem  are  negligible  for 
%  cable  X,  y  displacement  calculations. 

theta(l)  =  thetaS; 
tau(l)  =  tauS; 

else  if  Shape  ==  3||4;  %lnitial  tensions  provided  by  arbitrary  force  at  end  of  cable 
(substitution  for  ballast/float) 

%  Pode  Example  1:  TV(1)  =  40495  TH(1)  =  23140 
%  TV(1)  cannot  equal  zero  as  arctangent  would  be  undefined. 

TV(1)=  40495; 

TH(1)=  -23140; 

%  Because  of  flow  orientation,  the  arbitrary  horizontal  force  is 
%  positive  for  a  towed  cable.  The  horizontal  force  is  negative  for  a 
%  moored  cable.  When  the  final  results  are  graphed,  the  direction  of 
%  flow  will  be  the  same.  This  does  not  seem  to  have  an  impact  on  the 
%  results. 

tau(l)  =  sqrt(TV(l)^2  +  TH(1)^2); 
theta(l)  =  atan(TV(l)/TH(l)); 

x(l)  =  0; 
y(l)  =  0; 

%  Both  x(l)  and  y(l)  provide  normalization  for  the  array.  The  starting 
%  point  of  (0,0)  is  where  the  arbitrary  forces  act  on  the  end  of  the  cable, 
end 
end 
end 
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%  Calculate  the  eross-sectional  area  of  the  eable. 
AC  =  pi  *  DC  *  DC  /  4;  %m^2 


for  i  =  1  :n; 

%  Calculate  lift  and  drag  coefficients  of  the  eable. 

Cl(i)  =  Co  *  (sin(theta(i)))^2  *  cos(theta(i))  *  sign(theta(i)); 

Cd(i)  =  Co  *  (sin(theta(i)))^3  *  sign(theta(i))  +  Cf; 

%  Caleulate  the  new  tau  and  theta  coeffieients  for  each  successive  segment. 
df_tau(i)  =  ((g  *  AC  *  (DENC  -  DENW)  *  sin(theta(i)))  +  (1/2  *  DENW  *  (V^2)  *  DC  * 
(Cd(i)  *  cos(theta(i))  -  Cl(i)  *  sin(theta(i)))))  *  ds; 

tau(i+l)  =  tau(i)  +  df_tau(i); 

df_theta(i)  =  ((1  /  tau(i+l))  *  (g  *  AC  *  (DENC  -  DENW)  *  cos(theta(i))  -  1/2  *  DENW  * 
(V^2)  *  DC  *  (Cd(i)  *  sin(theta(i))  +  Cl(i)  *  cos(theta(i)))))  *  ds; 

theta(i+l)  =  theta(i)  +  df_theta(i); 

%  Caleulate  the  x  and  y  position  for  each  segment. 
x(i+l)  =  cos(theta(i))  *  ds  +  x(i); 

y(i+l)  =  sm(theta(i))  *  ds  +  y(i); 

%  Caleulate  the  new  horizontal  and  vertical  components  of  tension. 

TH(i)  =  tau(i)  *  cos(theta(i)); 

TV(i)  =  tau(i)  *  sin(theta(i)); 
end 

tauN  =  sqrt(TV.^2  +  TH.^2);  %  Total  eable  tension  in  N 
taulb  =  tauN  *  0.22482;  %  Total  eable  tension  in  Ibf 

thetarad  =  atan2(TV,TH);  %  Cable  angle  per  segment  in  radians. 

%  Corrections  so  data  is  properly  graphed  for  eaeh  case. 

%  Equivalency  is: 

%  Towed  eable  is  the  same  as  Shape  =  3  is  the  same  as  Density  =  2. 

%  Moored  cable  is  the  same  as  Shape  =  4  is  the  same  as  Density  =  1 . 

if  Shape  ==  3  %  (0,0)  is  where  cable  connects  to  buoy  on  seasurface. 
e  =  y(n+l); 

y  =  y-e; 

f  =  x(n+l); 
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X  =  x-f; 

thetadeg  =  rad2deg(thetarad)  +180; 
figure(l) 

plot(x,y,'rVLineWidth',  2, 'Marker', 'MarkerSize',  2); 
title('2-D  Position  of  Towed  Cable') 
xlabel  ('Horizontal  Displacement  (m)'); 
ylabel  ('Depth  Below  Seasurface  (m)'); 

else  if  Shape  ==  4  %  assumes  (0,0)  is  where  cable  connects  to  ground, 
e  =  y(n+l); 

y  =  y-e; 

f  =  x(n+l); 

X  =  x-f; 

thetadeg  =  rad2deg(thetarad)  +180; 
figure(l) 

plot(x,y,'r','LineWidth',  2, 'Marker', 'MarkerSize',  2); 
title('2-D  Position  of  Moored  Cable') 
xlabel  ('Horizontal  Displacement  (m)'); 
ylabel  ('Height  Above  Seafloor  (m)'); 

else  if  Density  ==  1  %  assumes  (0,0)  is  where  cable  connects  to  ground, 
e  =  y(n+l); 

y  =  y-e; 

f  =  x(n+l); 

X  =  x-f; 

thetadeg  =  rad2deg(thetarad)  +180; 
figure(l) 

plot(x,y,'r','LineWidth',  2, 'Marker', '.',  'MarkerSize',  2); 
title('2-D  Position  of  Moored  Cable') 
xlabel  ('Horizontal  Displacement  (m)'); 
ylabel  ('Height  Above  Seafloor  (m)'); 

else  if  Density  ==  2  %  (0,0)  is  where  cable  connects  to  buoy  on  seasurface. 
e  =  y(n+l); 

y  =  y-e; 

f  =  x(n+l); 

X  =  x-f; 

thetadeg  =  rad2deg(thetarad)  +180; 
figure(l) 

plot(x,y,'k',  'LineWidth',  2,'Marker', '.',  'MarkerSize',  2); 
hold  on 

title('2-D  Position  of  Towed  Cable') 
xlabel  ('Horizontal  Displacement  (m)'); 
ylabel  ('Depth  Below  Seasurface  (m)'); 
end 
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end 

end 

end 

%axis  square 
grid  on 
hold  on 

%  thetadeg  eorresponds  to  the  angle  relative  to  the  horizontal.  The  first 
%  angle  is  where  the  eable  eonneets  to  the  ballast/fioat. 

%  tauN  is  the  tension  in  the  eable.  The  first  tension  eorresponds  to  where 
%  the  eable  eonneets  to  the  ballast/float. 

%  Comment  out  lines  410-425  when  evaluating  multiple  eables. 

%  When  final  eable  is  being  evaluated,  uneomment  lines  410-425  to  add 
%  seafloor  or  seasurfaee.  Otherwise,  only  one  eable  graph  will  show  with 
%  the  seasurfaee  and  seafloor.  The  seasurfaee/seafloor  are  used  a  a 
%  referenee  point  to  show  the  user  that  the  eable  will  violate  a  boundary. 

%  MATLAB  does  not  ehange  its  ealeulations  if  the  eable  violates  a  boundary. 

z  =  -L:L/8;  %  Sets  limits  for  seafloor  and  seasurfaee  horizontal  boundaries. 

%  Default  set  to  1 .25  eable  seope  with  more  emphasis  on  the  left  side  of 
%  (0,0)  as  eurrent  flows  from  right  to  left. 


if  Shape  ==  3 

plot  (z,0,'b','LineWidth',  2,'Marker',  'MarkerSize',  6); 
else  if  Shape  ==  4 

plot  (z,0,'kVLineWidth',  2,'Marker', '.',  'MarkerSize',  2); 
else  if  Density  ==  1 

plot  (z,0,'k','LineWidth',  2,'Marker', '.',  'MarkerSize',  2); 
else  if  Density  ==  2 

plot  (z,0,'b','LineWidth',  2,'Marker',  'MarkerSize',  6); 
end 
end 
end 
end 

legend  Nitronie50B  SpeetralOOO  Seafloor  %ohange  legend  as  neeessary 
(seafloor/seasurfaee  is  always  last) 

%  legend  100m  10m  Im  Seafloor 

%  legend  17.46mm  15.88mm  14.29mm  12.70mm  Seafloor 
%  legend  Sphere  StreamlineBody  Seafloor 
%  legend  O.Okt  0.5kt  l.Okt  1.5kt  2.0kt  2.5kt  3.0kt  Seafloor 
%  legend  Nitronio5 OB  SpeetralOOO  Seafloor 
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APPENDIX  B.  SYNTACTIC  FOAM  SPECIFICATIONS 


High 

BZ  Grade  SuntsM 


Deep  Water  Foam  (Bathypelagic  zone) 


107  Frank  Mossberg  Drive 
Attleboro,  MA  02703 
TEL  (508)  226  -  3907  FAX  (608)  226  -3902 
tech@esvntactic.com  www.esyntactic.com 


Overview 

Syntactic  foams  are  a  unique  combination  of  hollow  spheres,  a  resin  matrix,  and  other  additives.  When 
combined  and  processed  properly,  the  constituents  form  a  lightweight  homogeneous  material  having  high 
compressive  and  hydrostatic  strengths.  This  yields  a  product  able  to  withstand  the  high  hydrostahc  pressures 

experienced  by  today’s  manned  and  unmanned  subsea  vehicles. 

The  BZ  Grade  of  material  is  specially  engineered,  high  performance 
syntactics  capable  of  operating  within  much  of  the  Bathypelagic  Zone 
(4,000  meters  or  1 3, 1 00  feet).  Due  to  the  nature  of  the  syntactic 
constituent  materials,  the  BZ  Grade  is  available  in  a  range  of  densities 
and  strengths.  The  product  is  cast  in  block,  sheet  or  rod  form  and  has 
isotropic  properties.  It  may  be  cut  or  machined  to  any  size  or  shape  to 
consistent  buoyancy  at  depth. 

Each  block  is  cast  as  a  single  unit,  but  may  be  cut  to  fit  the  application 
requirement.  Blocks  or  trimmed  parts  may  also  be  bonded  together  to 
form  a  larger  stmcture.  Standard  block  size  is  6  inch  x  12  inch  xl  2  inch, 
but  custom  sizes  and  shapes,  including  round  stock,  are  available. 

Properties 

Properties  provided  below  are  typical  for  the  cast  block: 


Product 

Density 

Ib/ft^ 

(g/cc) 

Service 

Pressure 

psi 

(Bar) 

Service 

Depth 

feet 

(meters) 

Compressive 

Strength 

psi 

(Mpa) 

Compressive 

Modulus 

ksi 

(Gpa) 

Hydrostatic 

Cmsh 

psi 

(Bar) 

Weight 

gain 

24  hours 
@  depth 

BZ-24 

24  ±2 

3,200 

7,250 

2,900 

170 

4,500 

3% 

(0.39  ±.03) 

(221) 

(2,200) 

(20.0) 

(1.17) 

(310) 

Max 

BZ-26 

26  ±2 

4,300 

9,750 

3,950 

180 

5,350 

3% 

(0.42  ±.03) 

(296) 

(2,980) 

(27.2) 

(1.24) 

(370) 

Max 

BZ-28 

28  ±2 

5,500 

12,500 

4,880 

210 

6,200 

3% 

(0.45  ±.03) 

(379) 

(3,810) 

(33.6) 

(1.45) 

(427) 

Max 

BZ-30 

30  ±2 

6,000 

13,600 

5,200 

240 

7,000 

3% 

(0.48  ±.03) 

(414) 

(4,150) 

(35.9) 

(1.65) 

(482) 

Max 
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APPENDIX  C.  SPECTRA  FIBER  SPECIFICATIONS 


Honeywell 


Honeywell  Spectra®  fiber  1000 
high-strength,  light-weight 
polyethylene  fiber 

Product  Information  Sheet 
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SPECTRA®  fiber  1000 

high-strength,  light-weight  polyethylene  fiber 


PRODUCT  DESCRIPTION:  Spectra  fiber  1000,  the  second  in  a  series  of  Spectra  fibers, 
was  developed  to  meet  customers’  needs  for  increased  performance.  It  is  available  in 
a  multitude  of  deniers  for  use  in  a  wide  range  of  applications. 


This  extended-chain  polyethylene  fiber  has  one  of  the  highest  strength-to-welght  ratios  of  any  man-made  fiber.  Spectra 
fiber  1000  has  a  tenacity  15  to  20  percent  higher  than  that  of  Spectra  fiber  900.  Spectra  fiber  is.  pound  for  pound,  15 
times  stronger  than  steel,  more  durable  than  polyester  and  has  a  specific  strength  that  is  up  to  40  percent  greater  than 
that  of  aramid  fiber  Specific  performance  is  dependent  upon  denier  and  filannent  count. 


Applicatiorrs: 

•  Police  and  military  ballistic  vests  and 
helmets 

•  Composite  armor  for  vehicles  and  aircraft 

•  Marine  lines  and  commercial  fishing  nets 
■  Industrial  cordage  and  slings 

•  Sailcloth,  fishing  line  and  other  sporting 
equipment 

•  Cut-protection  products 


Product  Benefits: 

•  Light  enough  to  float  (0.97  g/cc 
specific  gravity) 

■  High  resistance  to  chemicals, 
water  and  ultraviolet  light 

■  Excellent  vibration  damping 

•  Highly  resistant  to  flex  fatigue 


•  Low  coefficient  of  friction 

•  Good  resistance  to  abrasion 

•  Low  dielectric  constant  makes 
it  virtually  transparent  to  radar 


PHYSICAL  PROPERTIES 

(Nominal) 

Spectra'  fiber  1000 

WeightAJnit  Length  (Denier) 

215 

275 

375 

435 

650 

1300 

2600 

(Oecitex) 

239 

306 

417 

483 

722 

1444 

2888 

Ultimate  Tensile  Strength  (g/den) 

38 

36 

35 

34.5 

36 

35 

34 

(Gpa) 

3.25 

3.08 

3.00 

2.95 

3.08 

3.00 

2.91 

Breaking  Strength  (lbs.) 

18.0 

22.0 

29.0 

33.0 

52 

100 

195 

Modulus  (g/den) 

1320 

1320 

1200 

1180 

1175 

1150 

1135 

(Gpa) 

113 

113 

103 

101 

101 

98 

97 

Elongation  (%) 

29 

3.1 

3.1 

32 

3.5 

3.3 

3.5 

Density  (g/cc) 

0.97 

0,97 

097 

097 

0.97 

0,97 

0.97 

(Ibsiln^) 

0.035 

0.035 

0.035 

0.035 

0,035 

0.035 

0.035 

FlIamentAow 

60 

60 

60 

120 

120 

240 

480 

Filament  (dpO 

36 

4.6 

6.3 

3.6 

5.4 

5.4 

54 

For  more  information: 
Honeywell  Specialty  Materials 

101  Columbia  Road 
Morristown  NJ,  07962 
Toll  Free:  1-800-695-5969 
Europe: +32  16  39 12  19 
Asia:  +86  21  2894  2000 
WWW.  honeywell  com/spectra 
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INITIAL  DISTRIBUTION  LIST 


1.  Defense  Teehnieal  Information  Center 
Ft.  Belvoir,  Virginia 

2.  Dudley  Knox  Library 
Naval  Postgraduate  Sehool 
Monterey,  California 

3.  Professor  Fotis  Papoulias 
Naval  Postgraduate  Sehool 
Monterey,  California 

4.  Professor  Joseph  Riee 
Naval  Postgraduate  Sehool 
Monterey,  California 

5.  Professor  Andres  Larraza 
Naval  Postgraduate  Sehool 
Monterey,  California 

6.  Professor  Knox  Millsaps 
Naval  Postgraduate  Sehool 
Monterey,  California 

7.  Professor  John  Colosi 
Naval  Postgraduate  Sehool 
Monterey,  California 

8.  Professor  David  Ford 
Naval  Postgraduate  Sehool 
Monterey,  California 

9.  RADM  Gerald  Ellis  (Ret) 

Naval  Postgraduate  Sehool 
Monterey,  California 

10.  LT  Andrew  Hendrieksen 
Naval  Postgraduate  Sehool 
Monterey,  California 

1 1 .  CAPT  Mark  Kohlheim 
SSC  Paeifie 

San  Diego,  California 
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12.  Mark  Gillcrist 
SSC  Pacific 

San  Diego,  California 

13.  Bill  Marn 
SSC  Pacific 

San  Diego,  California 

14.  Robert  Creber 
SSC  Pacifie 

San  Diego,  California 

15.  Christopher  Fleteher 
SSC  Paeifie 

San  Diego,  California 

16.  Paul  Baxley 
SSC  Pacifie 

San  Diego,  California 

17.  CAPT  Miehael  Byman 
NUWC 

Newport,  Rhode  Island 

18.  Dana  Hesse 

Office  of  Naval  Researeh 
Washington,  DC 

19.  Thomas  Drake 
NSWC  Carderock 

West  Bethesda,  Maryland 

20.  Harry  Cox 
Loekheed  Martin 
Bethesda,  Maryland 

2 1 .  Edward  J ohnson 

Chief  of  Naval  Operations  (N871B) 
Washington,  DC 


58 


